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Abstract 

We calculate the electroweak corrections to the production of W + W~Z and 
ZZZ at the linear collider in the Standard Model. These processes are important 
for the extraction of the quartic couplings of the massive gauge bosons which can 
be a window on the mechanism of spontaneous symmetry breaking. We find that 
the weak corrections to some kinematic distributions show new features and hence 
cannot be explained by an overall scale factor. 
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1 Introduction 



Due to its clean environment an e + e~ linear collider in the TeV range is an ideal machine 
to probe in detail and with precision the inner working of the electroweak structure, 
in particular the mechansim of symmetry breaking. From this perspective the study of 
e + e~ — > W + W~Z and e + e~ — > ZZZ may be very instructive and would play a role 
similar to e + e~ — > W + W~ at lower energies. Indeed it has been stressed that e + e~ — > 
W + W~ Z and e + e~ — > ZZZ are prime processes for probing the quartic vector boson 
couplings [lj. In particular deviations from the gauge value in the quartic W + W~ZZ 
and ZZZZ couplings that are accessible in these reactions might be the residual effect 
of physics intimately related to electroweak symmetry breaking. Since these effects can 
be small and subtle, knowing these cross sections with high precision is mandatory. This 
calls for theoretical predictions taking into account loop corrections. 

Radiative corrections to e + e~ — > ZZZ have appeared recently in [2j and those to 
e + e~ — > W + W~Z in [3]. We have made an independent calculation of the electroweak 
corrections to e + e~ — > W + W~Z and e + e _ — > ZZZ, see [4|. Our preliminary results, 
eventually confirmed, have been presented in this workshop [5] prior to |3J. A detailed 
comparison between our results and the ones of Refs. [2J [3] has been done in [I]. 

In this report we summarize our results and make a further study on some distributions 
for W + W~Z production. 

2 Calculational details 

Our calculations are done in the framework of the SM. At leading order the process 
e + e~ — > ZZZ contains two types of couplings eeZ and ZZH. This process could probe 
the effect of a quartic ZZZZ coupling which is absent at tree-level, in the SM. The process 
e + e~ — > W + W~ Z is much more complicated with the involvement of trilinear and quartic 
gauge couplings in addition to the similar couplings as in ZZZ production. Compared 
to the well-tested process e + e~ — > W + W~ the new ingredients here are the two quartic 
gauge couplings WW Z^ and WW ZZ. Thus, this WW Z production at the ILC will be 
an excellent channel for studying these couplings. 

We have performed the calculation in at least two independent ways both for the 
virtual and the real corrections leading to two independent numerical codes (one code is 
written in Fortran 77, the other in C++). A comparison of both codes has shown full 
agreement at the level of the integrated cross sections as well as all the distributions that 
we have studied. 

Input parameters and renormalisation scheme: 

We follow closely the on-shell renormalisation scheme as detailed in Refs. [SI E]- To 
make the final results independent of the light quark masses we adopt a variant of 
the scheme. At tree level, the electromagnetic coupling constant is calculated as 
«g m = VSG^M^Sw/tc This absorbs some universal ml corrections and also the large 
logarithmic universal corrections proportional to ln(g 2 / m'j) where q is some typical en- 
ergy scale and mj a fermion mass. To avoid double counting we have to subtract the 
one-loop part of the universal correction from the explicit O(a) corrections by using the 
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counterterm 5Z e ^ = SZ e — Ar/2, the expression for Ar can be found in [7]. For one-loop 
corrections we use the coupling a(0) for both virtual and real photons. Thus the NLO 
corrections are of order a(0). Further details and the complete set of input parameters 
are given in |4]. 
Virtual corrections: 

The virtual corrections have been evaluated using a conventional Feynman-diagram based 
approach using standard techniques. We use the packages FeynArts and FormCalc-6 . to 
generate all Feynman diagrams and helicity amplitude expressions [8]. We also use SloopS 
to check the correctness of the amplitudes by checking non-linear gauge invariance (see 
[9] and references therein). The total number of diagrams in the 't Hooft-Feynman gauge 
is about 2700 including 109 pentagon diagrams for e + e~ — > W + W~Z and about 1800 
including 64 pentagons for e + e _ — > ZZZ. This already shows that e + e~ — > W + W~ Z 
with as many as 109 pentagons is more challenging than e + e~ — > ZZZ. Indeed get- 
ting stable results for all scalar and tensor (up to rank 4) box integrals in the process 
e + e~ — > W + W~ Z is a highly nontrivial task. An efficient way to solve this problem is 
by using higher precision arithmetic in part of the calculation. Further details related to 
loop integrals and the references for useful public codes are given in |4]. 
Real corrections: 

In addition to the virtual corrections we also have to consider real photon emission, i.e. 
the processes e + e~ — > W + W~ Z'y and e + e~ — > ZZZ'j. The corresponding amplitudes are 
divergent in the soft and collinear limits. The soft singularities cancel against the ones 
in the virtual corrections while the collinear singularities are regularized by the physical 
electron mass. To extract the singularities from the real corrections and combine them 
with the virtual contribution we apply both the dipole subtraction scheme and a phase 
space slicing method. The former is used to produce the final results since it yields smaller 
integration errors. Further details are given in |4]. 
Defining the weak corrections: 

It is well-known that the collinear QED correction related to initial state radiation in 
e + e~-processes is large. In order to see the effect of the weak corrections, one should 
separate this large QED correction from the full result. It means that we can define the 
weak correction as an infrared and collinear finite quantity. The definition we adopt in 
this paper is based on the dipole subtraction formalism. In this approach, the sum of the 
virtual and the so-called "endpoint" (see [TO] for the definition) contributions satisfies the 
above conditions and can be chosen as a definition for the weak correction 

'-'"weak C"virt ^endpoint • (1) 

For the numerical results shown in the next section, we will make use of this definition. 

3 Numerical results 

e+e~ ZZZ: 

As shown in Fig. [T] the tree-level cross section rises sharply once the threshold for pro- 
duction opens, reaches a peak of about 1.1 fb around a centre-of-mass energy of 600 GeV 
before very slowly decreasing with a value of about 0.9 fb at ITeV. The full NLO correc- 
tions are quite large and negative around threshold, —35%, decreasing sharply to stabilise 
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Figure 1: Left: the total cross section for e + e — > ZZZ as a function of a/s for the Born, 
full 0(a) and genuine weak correction. Right: the same for e + e~ — > W + W~Z . 



at a plateau around yfs = 600 GeV with —16% correction. The sharp rise and negative 
corrections at low energies are easily understood. They are essentially due to initial state 
radiation (ISR) and the behaviour of the tree-level cross section. The photon radiation 
reduces the effective centre-of-mass energy and therefore explains what is observed in the 
figure. On the other hand the genuine weak corrections, in the scheme, are relatively 
small at threshold, —7%. They however increase steadily with a correction as large as 
-18% at v^= ITeV. 
e+e- -> W+W-Z: 

Compared to ZZZ production, the cross section for e + e~ — > W + W~ Z is almost 2 orders 
of magnitudes larger for the same centre-of-mass energy. For example at 500 GeV it is 
about 40 fb at tree level, compared to 1 fb for the e + e~ — > ZZZ cross section. For an 
anticipated luminosity of lab -1 , this means that the cross section should be known at the 
per-mil level. The behaviour of the total cross section as a function of energy resembles 
that of e + e~ — > ZZZ. It rises sharply once the threshold for production opens, reaches a 
peak before very slowly decreasing as shown in Fig. [TJ However as already discussed the 
value of the peak is much larger, ~ 50 fb at NLO, moreover the peak is reached around 
yfs = 1 TeV, much higher than in ZZZ. This explains the bulk of the NLO corrections 
at lower energies which are dominated by the QED correction, large and negative around 
threshold and smaller at higher energies. As the energy increases the weak corrections 
get larger reaching about —18% at y/s = 1.5 TeV. 

In Fig. [2] we show the distributions in the rapidity and the transverse momentum of the 
W + . First, due to photon radiation, in the full NLO corrections some large corrections do 
show up at the edges of phase space. However when the QED corrections are subtracted, 
the weak corrections cannot be parameterized by an overall scale factor, for all the dis- 
tributions that we have studied. This feature together with possible normal/anomalous 
thresholds make it clear that calculating explicitly one-loop EW corrections is needed for 
a precise comparison with experimental data. 
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Figure 2: From top to bottom: distributions for the rapidity and the transverse momentum 
of the W + for e + e~ — > W + W~Z . The panels on the left show the tree-level, the full NLO 
and the weak correction. The panels on the right show the corresponding relative (to the 
tree-level) percentage corrections. 
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